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Introduction 


The  phosphoinositide-3-kinase  (PI3K)  pathway  regulates  numerous  critical  cellular 
functions  such  as  proliferation,  growth  and  survival  and  is  one  of  the  most  frequently 
altered  pathways  in  human  cancer  [1],  Tumor  cells  have  evolved  numerous  ways  of 
over-activating  this  pathway  to  induce  tumorigenesis,  such  as  amplifications  or 
mutations  in  upstream  receptor  tyrosine  kinases  (RTK),  mutations  or  deletion  of  the 
negative  regulator  PTEN,  somatic  mutations  in  PI3K  itself,  and  combinations  of  the 
above  [2-4],  Nearly  all  of  these  genetic  aberrations  have  been  detected  in  breast 
carcinomas.  This  supports  the  notion  that  activation  of  the  PI3K  and  the  nodal  kinase, 
protein  kinase  B  (AKT/PKB),  is  a  crucial  element  of  mammary  cell  viability.  However,  it 
also  implies  that  this  pathway  must  be  carefully  regulated  in  order  to  prevent 
uncontrolled  proliferation  and  survival  that  may  lead  to  breast  cancer. 

Identification  of  these  regulatory  modules  is  critical  to  the  development  of  novel 
treatments  and  to  better  inform  our  use  of  current  therapies  for  breast  cancers  driven  by 
PI3K  signaling.  We  have  now  uncovered  a  previously  uncharacterized  mechanism  to 
negatively  regulate  PI3K  signaling  in  breast  epithelial  cell  populations.  We  used  an 
approach  that  analyzes  signal  transduction  on  the  single  cell  level,  which  reveals 
population-context  dependent  regulatory  networks  that  cannot  be  detected  by  more 
traditional  analytical  means  such  as  western  blots.  Using  multi-color  flow  cytometry,  we 
demonstrate  that  MCF10A  mammary  cell  populations  invariably  maintain  heterogeneity 
in  PI3K  protein  levels,  characterized  by  a  bimodal  distribution  of  cells  expressing  either 
high  or  low  levels  of  the  enzyme.  We  further  show  that  modulation  of  PI3K  protein  levels 
is  regulated  by  proteosome-mediated  degradation  and  de  novo  synthesis.  This 
variability  in  PI3K  protein  levels  subsequently  leads  to  a  bimodal  distribution  in 
downstream  AKT  activity.  Our  data  lead  us  to  hypothesize  that  heterogeneous 
PI3K/AKT  activity  is  required  to  maintain  low  levels  of  PI3K  activity  in  large  cell 
populations  in  order  to  circumvent  oncogenesis  and  senescence  and  respond  to 
confluence  arrest. 

This  work  is  builds  on  a  growing  body  of  work  that  characterizes  cell-to-cell  variability  in 
cell  populations.  Ferrell  and  Machleder  demonstrated  that  variability  in  MAPK 
phosphorylation  leads  to  either  a  G2  or  metaphase  arrest  in  Xenopus  oocytes  [5]. 
Similarly,  high  or  low  levels  of  Nanog  determines  the  potential  of  embryonic  stem  cells  to 
terminally  differentiate  [6,  7],  and  high  or  low  levels  of  Sca-1  dictates  the  proclivity  for 
hematopoietic  progenitor  cells  to  commit  to  the  erythroid  or  the  myeloid  lineage  [8].  In 
certain  epithelial  cells,  gefitinib-resistance  can  be  conferred  by  high  IGF-IFt  signaling 
and  high  KDM5A  expression  [9]  and  camptothecin-resistance  by  high  DDX5  or  RFC1 
expression  [10].  These  studies  shed  much  light  on  the  importance  of  heterogeneous 
protein  expression  in  the  determination  of  distinct  cellular  fates.  Our  study  now  shows 
that  heterogeneity  in  signal  transduction  is  another  meaningful  source  of  cell-to-cell 
variability  that  results  not  in  distinct  cell  fates,  but  rather  influences  the  transient  behavior 
of  the  entire  population.  This  is  particularly  important  in  the  context  of  tumor  growth  and 
the  maintenance  of  homeostasis  in  tissues  and  organs. 
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Body 


We  present  this  body  of  work  in  its  entirety  in  the  Appendix  A.  We  are  pleased  to  report 
that  it  has  recently  been  accepted  for  publication  in  Current  Biology  (Cell  Press).  The 
major  findings  are  outlined  briefly  below.  The  first  five  findings  were  described  in  the 
previous  Annual  Summary,  and  the  remaining  three  were  completed  in  the  time  since 
the  submission  of  the  last  Annual  Summary. 

AKT  activation  in  MCF10A  populations  is  heterogeneous  and  exhibits  a  bimodal 
distribution 

Upon  acute  stimulation  with  EGF,  we  observed  striking  heterogeneity  in  AKT  activation 
characterized  by  the  bimodal  distribution  of  cells  into  a  pAKT-negative  population 
(-70%)  or  a  pAKT-positive  population  (-30%).  Bimodality  was  also  observed  following 
insulin  stimulation,  suggesting  that  this  is  a  general,  rather  than  ligand-specific,  effect  on 
PI3K  activation.  We  also  observed  robust  bimodal  activation  of  AKT  in  another  non- 
transformed  human  mammary  epithelial  cell  line,  HMECs,  but  not  in  fully  transformed 
breast  cancer  cell  lines.  This  suggests  that  cancer  cells  have  evolved  a  mechanism  to 
overcome  this  negative  regulation  of  AKT  activation. 

pi  10a  protein  level  is  variable  in  single  cells  and  determines  AKT  activity 

We  found  that  pi  10a  protein  levels,  like  pAKT,  exhibit  bimodality  in  cell  populations, 
pi  10a  protein  level  also  positively  correlated  with  pAKT  status  (R2  =  0.68),  indicating 
that  only  cells  with  high  levels  of  pi  10a  can  activate  AKT.  Additionally,  levels  of  p85a 
mirror  that  of  pi  10a  in  single  cells,  confirming  that  the  entire  holoenzyme  is  required  to 
induce  downstream  signaling. 

pi  10a  protein  levels  are  dynamic  within  single  cells 

We  observed  dynamic  changes  in  pi  10a  protein  levels  following  EGF  stimulation,  which 
had  not  been  previously  reported.  We  show  that  pi  10a  protein  in  pAKT-positive  cells 
undergoes  a  precipitous  drop  10  minutes  after  EGF  stimulation,  when  AKT  activity  is 
maximal.  In  contrast,  the  major  subgroup  of  cells  with  low  levels  of  pAKT  does  not  show 
a  significant  decline  in  pi  10a  during  the  time  course  of  EGF  stimulation. 

Dynamic  changes  in  pi  10a  protein  levels  are  regulated  by  a  degradation/re¬ 
synthesis  cycle 

We  treated  cells  with  the  proteosome  inhibitor,  MG132,  or  the  translation  inhibitor, 
cyclohexamide,  to  monitor  effects  on  pi  10a  protein  levels.  We  demonstrate  that 
proteosome-mediated  degradation  of  pi  10a  transitions  cells  from  the  pi  1  oahl9h/medlum 
state  to  the  pi  10alow  state,  thereby  negatively  regulating  AKT  activity.  Additionally,  we 
show  that  pi  1 0a  re-synthesis  is  responsible  for  shifting  cells  from  the  pi  1 0alow  state 
back  to  pi  10ahigh  state,  thereby  generating  cells  competent  to  activate  AKT. 

p110ah'9h  cells  initiate  colony  formation 

We  seeded  single  cells  to  generate  clonal  populations  and  show  that  early  passage 
clones  exhibit  the  reversed  bimodal  profile,  characterized  by  a  major  proportion  of  cells 
with  high  levels  of  pi  10/pAKT  and  a  small  proportion  of  cells  with  low  levels  of 
pi  1 0/pAKT.  These  results  suggest  that  only  pi  1 0ahigh  cells  were  capable  of  forming 
colonies.  Interestingly,  as  these  clones  were  serially  passaged,  they  all  reverted  back  to 
the  parental  bimodal  distribution.  This  suggests  that  pi  1 0ah'9h  cells  are  critical  for  colony 
formation  and  early  clonal  expansion,  however  heterogeneous  populations  of 
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predominantly  pi  1 0alow  cells  facilitate  exponential  growth.  Populations  that  fail  to  revert 
undergo  senescence. 

Hotspot  PIK3CA  mutations  confer  stabilization  of  the  p110ah'9h  state 

We  analyzed  MCF10A  cells  expressing  wildtype  or  either  of  two  common  oncogenic 
pi  10a  mutants,  H1047R  and  E545K.  We  show  that  the  mutant  proteins  are  maintained 
at  higher  concentrations  in  cells.  This  stability  is  partially  conferred  by  a  resistance  to 
degradation  and  results  in  a  greater  proportion  of  cells  with  high  levels  of  pi  10a,  which 
ultimately  leads  to  much  higher  levels  of  AKT  activity  within  the  population. 

Establishing  cell-cell  contacts  stabilizes  the  pi  10alow  state 

We  measured  pi  10a  levels  of  exponentially  growing  populations  every  24h  and  found 
that  pi  10a  and  pAKT  levels  steadily  decreased  as  cell  density  increased.  Under  low- 
density  growth  conditions,  the  pi  1 0ahigh  population  is  enriched  for  cells  with  0-1  cell-cell 
contacts,  while  the  pi  1 0alow  population  is  enriched  for  cells  with  2-4  cell-cell  contacts. 
Thus,  cell-cell  contacts  may  modulate  pi  10a  levels  in  order  to  establish  contact 
inhibition/confluence  arrest. 

AKT  activity  is  enriched  at  the  unencumbered  edges  of  tumors 

We  stained  mammary  tumors  arising  from  MMTV-Cre;  BRCA1  flox/flox;  p53+/-  mice  for 
pAKT  S473.  We  detected  heterogeneous  pAKT-positivity  throughout  the  tumor, 
including  an  enrichment  of  pAKT-positive  cells  along  pushing  margins  and  in  the  lumens 
of  ductal  hyperplasias.  Both  of  these  regions  represent  areas  with  fewer  cell-cell 
contacts,  supporting  our  hypothesis  that  low  cell  density  enhances  the  pAKT-positive 
state. 
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Key  Research  Accomplishments 

We  have  conducted  an  investigation  of  PI3K  signaling  in  normal  and  cancer  breast 
epithelial  cell  populations  using  multi-color  flow  cytometry.  We  have  characterized  a 
novel  mechanism  to  negatively  regulate  PI3K  signaling  in  cell  populations  that  is  based 
on  the  maintenance  of  variable  PI3K  protein  expression.  The  key  findings  are: 

•  AKT  activation  in  MCF10A  populations  is  heterogeneous  and  exhibits  a  bimodal 
distribution 

•  pi  10a  protein  level  is  variable  in  single  cells  and  determines  AKT  activity 

•  pi  10a  protein  levels  are  dynamic  within  single  cells 

•  Dynamic  changes  in  pi  1 0a  protein  levels  are  regulated  by  a  degradation/re¬ 
synthesis  cycle 

•  pi  1 0ah'9h  cells  initiate  colony  formation 

•  Hotspot  PIK3CA  mutations  confer  stabilization  of  the  pi  1 0ah'9h  state 

•  Establishing  cell-cell  contacts  stabilizes  the  pi  1 0alow state 

•  AKT  activity  is  enriched  at  the  unencumbered  edges  of  tumors 
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Reportable  Outcomes 


During  the  funding  period  of  this  BCRP  award,  the  following  manuscripts  have  been 
published  or  have  been  accepted  for  publication: 

1)  Yuan,  T.L.,  and  Cantley,  L.C.  Introduction.  CurrTop  Microbiol  Immunol  346,  1-7. 

2)  Yuan,  T.L.,  Wulf,  G.,  Burga,  L.,  and  Cantley,  L.C.  (2011)  Cell-to-cell  variability  in 
PI3K  protein  regulates  PI3K-AKT  pathway  activity  in  cell  populations.  Curr  Biol. 
Accepted. 

This  work  was  presented  in  poster  format  at  the  2010  Cold  Spring  Harbor  James 
Watson  Cancer  Symposium. 

This  work  has  been  presented  in  seminar  format  at  Genentech  and  UCSF. 
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Conclusions 


Single  cell  analysis  of  signal  transduction  networks  can  uncover  many  phenomena  in 
cell  populations  that  techniques  that  measure  population  averages  cannot.  Cell-to-cell 
heterogeneity  is  one  such  phenomenon.  Heterogeneity  can  be  the  result  of  genetic  or 
non-genetic  variation  and  can  have  major  effects  on  the  behavior  of  the  entire 
population.  We  have  shown  that  heterogeneity  in  PI3K  protein  and  pathway  activation  is 
an  important  mechanism  in  mammary  epithelial  cells  to  regulate  growth  in  variable 
environmental  contexts.  Under  acute  growth  factor  stimulation,  we  show  that 
heterogeneity  results  in  AKT  activation  in  only  a  subset  of  cells,  which  ensures  that  PI3K 
activity  does  not  reach  oncogenic  levels.  In  low  cell  density  conditions,  we  show  that 
cells  with  high  levels  of  PI3K  are  selected  for  colony  formation  and  outgrowth,  while 
under  high  density  conditions,  cells  are  maintained  at  the  low  PI3K  state  to  dampen 
growth  signals  to  perhaps  to  induce  confluence  arrest. 

In  cancer  cells  this  negative  regulation  appears  to  be  overcome.  Cells  expressing 
mutant  forms  of  pi  10  appear  to  do  this  by  stabilizing  the  PI3K  enzyme,  thus  maintaining 
more  actively  signaling  cells  within  the  population.  However,  we  showed  that  breast 
carcinomas  derived  from  MMTV-Cre;BRCA1flox/flox;p53+/'  mice  still  maintain  heterogeneity 
in  AKT  activation.  There  are  likely  to  be  many  more  AKT-positive  cells  than  in  normal 
breast  tissue,  however  this  remaining  heterogeneity  posits  a  major  obstacle  to  the 
efficacy  of  PI3K  inhibitors  in  the  treatment  of  breast  cancer. 

It  is  likely  that  only  tumor  cells  with  high  levels  of  pathway  activation  will  die  in  response 
to  acute  inhibition  of  PI3K.  Therefore,  heterogeneous  tumors  will  exhibit  only  “fractional 
killing”  in  response  to  a  PI3K  or  AKT  inhibitor.  However,  it  may  be  possible  to  kill  a 
larger  fraction  of  the  tumor  cells  by  adjusting  the  therapy  to  ensure  that  the  subset  of 
cells  that  are  protected  from  PI3K  inhibitor  therapy  at  the  time  of  the  initial  dose,  become 
exposed  to  the  inhibitor  a  second  time  as  they  cycle  into  a  high  PI3K  pathway  state. 
Thus,  a  detailed  understanding  of  the  pharmacokinetics  of  the  drug  is  important  in 
deciding  the  intensity  and  frequency  of  dosing.  While  continuous  high  doses  of  PI3K 
inhibitors  over  several  days  or  weeks  might  effectively  kill  all  tumor  cells  as  they  cycle 
into  the  high  activity  state,  dose-limiting  toxicities  might  preclude  such  treatment. 

A  number  of  pre-clinical  studies  using  PI3K  inhibitors  may  support  this  hypothesis.  Once 
daily  oral  administration  of  a  PI3K  inhibitor  to  tumor-bearing  mice  results  in  delayed 
tumor  growth  but  rarely  tumor  regression.  pAKT  is  ablated  in  these  tumors  1-2h  after 
dosing,  but  reemerges  several  hours  later  [1 1  -1 5].  Our  data  would  suggest  that  cells 
that  were  PI3Klow;pAKTlow  during  the  initial  dose  escaped  the  toxic  effects  of  the  drug, 
and  cycled  to  the  PI3Khigh;pAKTh'9h  state  several  hours  later  when  concentrations  of  the 
drug  decreased.  The  possibility  that  PI3K  inhibitors  are  selectively  targeting  the 
advancing  edges  or  “pushing  margins”  of  tumors  is  also  consistent  with  fact  that  these 
tumors  are  not  advancing,  but  rather  stabilizing. 

We  believe  that  this  work  brings  to  light  the  need  to  carefully  match  the  degree  of  tumor 
heterogeneity  with  the  appropriate  dosage  of  a  drug  in  order  to  achieve  sufficient  tumor 
regression. 
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Summary 

Background 

Cell-to-cell  variability  in  populations  has  been  widely  observed  in  mammalian  cells.  This 
heterogeneity  can  result  from  random  stochastic  events  or  can  be  deliberately 
maintained  through  regulatory  processes.  In  the  latter  case,  heterogeneity  should  confer 
a  selective  advantage  that  benefits  the  entire  population. 

Results 

Using  multicolor  flow  cytometry,  we  have  uncovered  robust  heterogeneity  in  PI3K  activity 
in  MCF10A  cell  populations,  which  had  been  previously  masked  by  techniques  that  only 
measure  population  averages.  We  show  that  AKT  activity  is  bimodal  in  response  to 
EGF  stimulation  and  correlates  with  PI3K  protein  level,  such  that  only  cells  with  high 
PI3K  protein  can  activate  AKT.  We  further  show  that  heterogeneity  in  PI3K  protein 
levels  is  invariably  maintained  in  cell  populations  through  a  degradation/re-synthesis 
cycle  that  can  be  regulated  by  cell  density. 

Conclusions 

Given  that  the  PI3K  pathway  is  one  of  the  most  frequently  altered  pathways  in  cancer, 
we  propose  that  heterogeneity  in  PI3K  activity  is  beneficial  to  normal  tissues  by 
restricting  PI3K  activation  to  only  a  subset  of  cells.  This  may  serve  to  protect  the 
population  as  a  whole  from  over-activating  the  pathway,  which  can  lead  to  cellular 
senescence  or  cancer.  Consistent  with  this,  we  show  that  oncogenic  mutations  in  pi  10a 
(H1047R  and  E545K)  partially  evade  this  negative  regulation,  resulting  in  increased  AKT 
activity  in  the  population. 

Highlights 

•  AKT  activity  in  MCF10A  populations  is  heterogeneous  and  is  regulated  by  PI3K 
protein  level 

•  Dynamic  changes  in  PI3K  protein  levels  are  regulated  by  degradation  and  de 
novo  translation 

•  Cells  with  high  and  low  levels  of  PI3K  protein  have  distinct  cellular  functions 

•  Oncogenic  PIK3CA  mutations  stabilize  high  PI3K  expression  and  induce  high 
AKT  activity 
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Introduction 


Among  the  thousands  of  proteins  expressed  in  a  single  cell  at  any  one  time,  it  is 
inevitable  that  even  genetically  identical  cells  will  have  variable  expression  of  certain 
proteins.  This  cell-to-cell  heterogeneity  can  have  profound  effects  on  the  phenotype  of  a 
population.  In  bacteria  and  yeast,  a  heterogeneous  population  has  a  better  chance  of 
surviving  in  fluctuating  environments  [1-3].  In  cancer  cells,  variable  protein  expression 
results  in  non-uniform  responses  to  anti-cancer  agents  and  only  “fractional”  killing,  which 
is  beneficial  to  the  tumor  but  detrimental  to  the  patient  [4-7], 

Population  heterogeneity  can  be  the  result  of  genetic  and  non-genetic  variability  [8]. 
Non-genetic  sources  of  heterogeneity  are  of  particular  interest  because  they  are  not 
fixed  in  the  population  by  heritable  transmission.  Rather,  these  sources  can  be 
stochastic  in  nature,  such  as  random  fluctuations  in  transcription  or  translation,  or  they 
can  be  regulated,  such  as  through  feedback  mechanisms,  so  that  variability  is 
deliberately  maintained.  Flow  cytometry  allows  us  to  measure  protein  concentrations  in 
single  cells  within  large  populations,  which  has  proved  useful  in  differentiating  between 
stochastic  and  regulated  sources  of  variation.  Variable  protein  levels  caused  by 
stochastic  events  are  usually  rare  and  slightly  increase  the  coefficient  of  variance  (CV)  of 
an  otherwise  log-normal  concentration  distribution  [9].  Regulated  variance  usually 
deviates  from  tight  log-normal  distributions  and  in  the  most  extreme  cases  might  display 
bimodality  [10].  Additionally,  variance  in  the  expression  level  of  these  proteins  tend  to 
correlate  with  variance  in  functionally  related  proteins  [9]. 

We  demonstrate  cell-to-cell  variability  in  phosphoinositide  3-kinase  (PI3K)  pathway 
activation  in  mammary  epithelial  cells  that  is  indeed  both  bimodal  and  correlates  with 
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multiple  proteins  in  the  pathway.  This  implies  that  heterogeneous  PI3K  activation  in  cell 
populations  is  maintained  by  a  regulated  process  and  likely  confers  a  selective 
advantage.  The  PI3K  pathway  regulates  numerous  critical  cellular  functions  such  as 
proliferation,  growth  and  survival  and  is  one  of  the  most  frequently  altered  pathways  in 
human  cancer  [11].  Tumor  cells  have  evolved  numerous  ways  of  over-activating  this 
pathway  to  induce  tumorigenesis,  such  as  amplifications  or  mutations  in  upstream 
receptor  tyrosine  kinases  (RTK),  mutations  or  deletion  of  the  negative  regulator  PTEN, 
somatic  mutations  in  PI3K  itself,  and  combinations  of  the  above  [12-14].  In  this  view,  we 
propose  that  cell-to-cell  variability  in  PI3K  activity  serves  as  a  mechanism  to  keep  overall 
PI3K  activity  limited  in  normal  tissues  to  avoid  hyperactivity  that  may  lead  to  cancer. 
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Results 


AKT  activation  in  MCF10A  populations  is  heterogeneous  and  exhibits  a  bimodal 
distribution 

To  assess  cell-to-cell  variability  in  PI3K  pathway  activation,  serum  and  growth  factor 
starved  MCF10A  cells  were  acutely  stimulated  with  EGF,  and  AKT  activity  was 
measured  in  single  cells  by  flow  cytometry.  As  expected,  in  the  absence  of  EGF  or  with 
pretreatment  with  the  PI3K  inhibitor  wortmannin,  cells  exhibited  a  log-normal  distribution 
in  the  pAKT-negative  gate  (Figure  1  A).  Following  treatment  with  saturating  amounts  of 
EGF,  we  observed  striking  heterogeneity  in  AKT  activation  characterized  by  the  bimodal 
distribution  of  cells  into  the  pAKT-negative  (-70%)  or  pAKT-positive  gates  (-30%) 

(Figure  1A).  The  percentage  of  cells  in  the  positive  gate  was  maximal  at  5-10  minutes 
and  decreased  over  time.  We  found  that  the  bimodality  was  not  due  to  differential  EGFR 
activation  or  total  AKT  protein  abundance,  both  of  which  exhibit  log-normal  distributions 
and  show  weak  correlation  with  pAKT  status  (Figure  SI  A,B).  Bimodality  was  also 
observed  following  insulin  stimulation,  suggesting  that  this  is  a  general,  rather  than 
ligand-specific,  effect  on  PI3K  activation  (Figure  SIC).  We  also  observed  robust 
bimodal  activation  of  AKT  in  another  non-transformed  human  mammary  epithelial  cell 
line,  PIMECs  (Figure  S1E). 

As  reported  in  the  literature,  cell  populations  expressing  the  oncogenic  FI1047R  PIK3CA 
mutation  have  higher  AKT  activity  [1 5-1 7].  The  enzyme  encoded  by  this  mutant  form  of 
PIK3CA  has  been  shown  to  have  higher  specific  activity  than  that  encoded  by  wild-type 
PIK3CA  [15,  16].  To  explore  the  possibility  that  this  mutation  shifts  cells  uniformly  to  the 
pAKT-positive  state,  MCF10A  cells  stably  expressing  this  mutant  form  of  pi  10a  were 
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analyzed  by  flow  cytometry.  We  observed  that  even  in  the  presence  of  oncogenic 
pi  10a,  cells  still  segregate  into  a  bimodal  distribution  and  retain  a  substantial  population 
of  non-responders.  The  H1047R  mutation  in  MCF10A  cells  thus  achieves  higher 
average  pAKT  levels  by  shifting  more  cells  into  the  pAKT-positive  state  (Figure  1  A). 

This  implies  that  the  mechanism  that  maintains  heterogeneity  in  cell  populations  is  both 
robust  and  is  not  eliminated  by  expression  of  a  mutant  form  of  PI3K  with  constitutively 
high  specific  activity. 

Lastly,  we  wanted  to  ensure  that  heterogeneous  AKT  activation  is  preserved  in  more 
physiological  settings.  Debnath  et  al  have  previously  observed  sporadic  pAKT  staining 
in  MCF10A  3D  acinar  structures,  which  is  important  for  proper  lumen  formation  [18,  19]. 
To  see  if  this  heterogeneous  AKT  activation  was  observed  in  vivo ,  we  analyzed 
spontaneously  arising  tumors  from  MMTV-Cre;  BRCA1  flox/flox;  p53+/-  mice  and  indeed 
observed  heterogeneity  in  pAKT  positivity  (Figure  IB). 

pi  10a  protein  level  is  variable  in  single  cells  and  determines  AKT  activity 

Given  that  the  bimodal  AKT  activation  could  not  be  explained  by  bimodal  distribution  of 
activated  EGFR  or  total  AKT  protein,  we  next  looked  for  variability  in  PI3K  protein 
abundance.  There  are  no  pi  10a  antibodies  amenable  for  flow  cytometry,  however 
pAKT  bimodality  was  maintained  in  MCF1 0A  cells  stably  expressing  FIA-tagged  pi  1 0a. 
We  thus  utilized  anti-FIA  antibodies  to  monitor  levels  of  pi  10a.  Using  single  cell 
analysis,  we  found  that  pi  10a  protein  levels,  like  pAKT,  exhibit  bimodality  in  cell 
populations  (Figure  1C),  pi  10a  protein  level  also  positively  correlated  with  pAKT  status 
(R2  =  0.68),  indicating  that  only  cells  with  high  levels  of  pi  10a  can  activate  AKT.  Cells 
with  medium  levels  of  pi  10a  have  intermediate  levels  of  AKT  activity,  suggesting  that 
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these  cells  may  be  transitioning  dynamically  between  the  two  dominant  modes  (Figure 
1C).  Subpopulations  in  these  varying  states  can  also  be  visualized  by 
immunofluorescence  (Figure  ID). 

The  p85a-p1 10a  heterodimer  is  a  very  stable  complex  that  does  not  readily  dissociate, 
and  it  is  thus  thought  that  pi  1 0a  levels  generally  correlate  with  p85a  levels  [20,  21  ].  We 
monitored  the  distribution  of  endogenous  p85a  in  MCF10A  cells  and  show  that  levels  of 
endogenous  p85a  mirror  that  of  FIA-tagged  pi  1 0a  in  single  cells  (Figure  SID).  This 
correlation  indicates  that  the  heterogeneity  we  observe  is  not  an  artifact  of  exogenously 
expressed  pi  10a  and  confirms  that  high  levels  of  the  entire  functional  holoenzyme  are 
required  for  AKT  activation,  as  expected. 

pi  10a  protein  levels  are  dynamic  within  single  cells 

In  agreement  with  previous  studies  [22-24],  total  PI3K  levels  do  not  appear  to  change 
upon  growth  factor  stimulation  as  measured  by  western  blot  (Figure  2A).  To  test  the 
possibility  that  pi  1 0a  protein  levels  are  dynamic  on  the  subpopulation  level,  we 
measured  pi  10a  levels  in  pAKT-positive  and  pAKT-negative  subpopulations  over  a  time 
course  of  EGF  stimulation  by  single  cell  analysis.  Surprisingly,  we  observed  dynamic 
changes  in  pi  10a  protein  in  pAKT-positive  cells.  This  minor  subgroup  of  cells  has  the 
highest  levels  of  pi  10a  (and  highest  levels  of  pAKT)  and  undergoes  a  precipitous  drop 
in  pi  10a  protein  at  10  minutes  after  EGF  stimulation,  when  AKT  activity  is  maximal 
(Figure  2B).  Over  the  next  4  hours,  this  population  gradually  recovers  pi  10a  levels 
(Figure  2B).  In  contrast,  the  major  subgroup  of  cells  with  low  levels  of  pAKT  does  not 
show  a  significant  decline  in  pi  1 0a  (and  do  not  show  activation  of  AKT)  during  the  time 
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course  of  EGF  stimulation.  We  attempted  to  verify  these  results  with  live  cell  imaging  of 
fluorescently  tagged  pi  1 Oa,  however  the  size  of  the  fluorescent  protein-pi  1 0  fusion 
prohibited  efficient  transfection  and  expression. 

However,  our  observed  changes  in  pi  10a  levels  also  correlate  with  a  striking  change  in 
the  proportion  of  cells  in  the  pi  1 0ahl9h,  pi  1 0amedlum,  and  pi  1 0alow  subpopulations  over 
the  time  course  of  EGF  stimulation.  After  10  minutes  of  EGF  stimulation,  there  is  a 
major  drop  in  the  proportion  of  cells  in  the  pi  10ah'9h  state  and  an  increase  in  the 
proportion  of  cells  in  the  pi  10amedlum  state  (Figure  2C).  Over  the  next  4  hours,  the 
populations  recover  to  the  original  distribution  (Figure  2C).  These  data  indicate  that 
there  is  a  significant  decrease  in  pi  1 Oct  levels  in  a  subpopulation  of  cells  characterized 
by  high  pAKT  and  high  pi  1 0ct.  The  failure  to  detect  this  acute  change  in  pi  1 0ct  by 
western  blot  can  be  explained  by  the  fact  that  this  population  represents  only  a  minor 
fraction  of  total  cells,  whereas  the  protein  in  a  western  blot  is  dominated  by  a  major 
fraction  of  cells  (70  to  80%)  that  fail  to  activate  the  PI3K/AKT  pathway  and  maintain 
static  pi  1 0ct  levels.  Endogenous  p85ct  levels  exhibit  nearly  identical  kinetics  to 
pi  1 0ct,  indicating  that  the  entire  holoenzyme  is  dynamically  regulated  to  modulate  AKT 
activity  (Figure  S2). 

We  also  plotted  the  time  course  of  EGF-stimulated  AKT  phosphorylation  in  the  sub¬ 
populations  of  cells  that  had  high,  medium  or  low  levels  of  pi  1 0ct  (Figure  2D).  The 
subset  of  cells  with  the  highest  level  of  pi  1 0a  not  only  had  the  highest  pAKT  signal 
amplitude,  but  also  the  longest  pAKT  signal  duration.  However,  the  overall  population 
response  was  dominated  by  the  major  subset  of  cells  that  had  minimal  AKT 
phosphorylation. 
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Dynamic  changes  in  pi  10a  protein  levels  are  regulated  by  a  degradation/re¬ 
synthesis  cycle 

We  next  sought  to  understand  the  mechanism  by  which  pi  1 0a  levels  were  modulated. 

In  these  cells,  HA-pl  1 0a  expression  is  driven  by  the  CMV  promoter,  suggesting  that 
variability  in  pi  10a  expression  is  post-transcriptionally  regulated.  To  test  if  the  drop  in 
pi  10a  levels  in  actively  signaling  cells  was  due  to  proteosome-mediated  degradation, 
we  treated  cells  with  MG  132  and  observed  an  increase  in  pi  10a  protein  on  the 
population  level  (Figure  3A).  On  the  subpopulation  level,  MG132  treatment  resulted  in 
the  accumulation  of  cells  in  the  pi  10amedlum  state  and  a  loss  of  cells  from  the  pi  10alow 
state  (Figure  3B).  This  suggests  that  proteosome-mediated  degradation  of  pi  10a 
transitions  cells  from  the  pi  1 0ahigh/medlum  state  to  the  pi  10alow  state,  thereby  negatively 
regulating  AKT  activity.  Short-term  treatment  of  cells  with  MG132  prior  to  acute  EGF 
stimulation  partially  rescues  the  drop  in  pi  1 0  levels  in  pAKT  high  cells  and  prevents  the 
shift  of  cells  from  the  pi  1 0ahigh  to  pi  1 0amedlum  state  (Figure  S3A,B). 

To  show  that  pi  10a  can  be  resynthesized  and  thus  transition  cells  from  the  pi  10alow  to 
the  pi  1 0amedlum/high  state,  we  treated  cells  with  cyclohexamide  (CFIX).  On  the  population 
level,  we  observed  an  overall  decrease  in  pi  10a  protein  with  CFIX  treatment  (Figure  3A). 
On  the  subpopulation  level,  we  observed  an  accumulation  of  cells  in  the  pi  10alow  state 
with  concomitant  loss  of  cells  from  the  pi  1 0amedlum  and  pi  1 0ahigh  states  (Figure  3B). 
These  data  show  that  pi  1 0a  re-synthesis  is  responsible  for  shifting  cells  from  the 
pi  1 0alow  state  to  pi  1 0ahigh  state,  thereby  generating  cells  competent  to  activate  AKT. 
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Our  results  thus  far  suggest  a  model  whereby  pi  10a  levels  in  single  cells  oscillate  via  a 
degradation/re-synthesis  cycle,  which  transitions  cells  through  phases  of  high  pi  10a 
expression  (competent  to  activate  the  pathway)  and  low  pi  10a  (incompetent  to  activate). 
Large  populations  contain  cells  in  various  stages  of  this  cycle,  thus  accounting  for  the 
vast  heterogeneity  we  observe  (Figure  3C).  Though  many  extrinsic  factors  are  likely  to 
influence  the  duration  and  progression  of  this  cycle,  we  found  no  correlation  between 
pi  10a  or  pAKT  status  with  cell  cycle,  as  assessed  by  DNA  content,  or  with  cell  size  or 
granularity  (Table  SI  and  Figure  S3C). 

Hotspot  PIK3CA  mutations  confer  stabilization  of  the  p110ah'9h  state 

Escaping  this  negative  regulation  of  AKT  may  be  particularly  important  in  cancer  cells, 
given  that  FI1047R-expressing  MCF10A  cells  and  several  tumor-derived  breast  cancer 
cell  lines  can  activate  AKT  in  a  greater  proportion  of  cells  compared  to  non-transformed 
cell  lines  (Figures  1A,  S1E).  To  test  if  this  was  due  to  the  ability  to  stabilize  the  pi  10ah'9h 
state,  we  measured  the  level  of  pi  10a  in  MCF10A  cells  stably  expressing  exogenous 
FIA-tagged  wildtype,  FI1047R  mutant,  or  E545K  mutant  pi  10a. 

Both  mutant  forms  of  pi  10a  were  expressed  at  higher  levels  than  wildtype  pi  10a 
(Figure  4A),  which  was  due  to  the  maintenance  of  more  cells  in  the  pi  10ah'9h  state  and 
fewer  cells  in  the  pi  10alow  state  (Figure  4B,C,H).  To  ensure  that  these  differences  were 
not  due  to  variability  in  retroviral  infection  efficiency,  multiple  batches  of  cell  lines  were 
tested,  and  all  lines  were  grown  under  selection  for  several  generations  before  use  in 
experiments. 
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We  next  sought  to  determine  if  the  mutant  forms  of  pi  1 0a  were  stabilizing  the  pi  1 0ah'9h 
state  by  resisting  degradation.  We  monitored  the  degradation  kinetics  of  the  mutant  cell 
lines  as  we  did  in  Figure  2B.  Though  both  mutant  cell  lines  underwent  rapid  pi  1 0a 
degradation  following  EGF  stimulation,  the  E545K  mutant  displayed  more  muted 
degradation  (though  this  was,  in  part  due  to  higher  basal  expression)  and  both  mutants 
exhibited  considerable  variability  between  experiments  compared  to  wildtype  (Figure 
4D).  Additionally,  MG132  treatment  of  asynchronously  growing  mutant  cells  did  not 
result  in  an  increase  in  overall  pi  10a  levels  as  observed  in  the  wildtype  line  (Figure  4E). 
These  results  suggest  that  the  FI1047R  and  E545K  mutations  confer  resistance  to 
proteosome-mediated  degradation  under  exponential  growth  conditions  that  may  be 
distinct  from  the  acute  degradation  that  occurs  following  EGF  stimulation  of  starved 
cells.  Notably,  all  cell  lines  responded  similarly  upon  CFIX  treatment,  ruling  out  the 
possibility  that  mutant  pi  10a  is  preferentially  synthesized  (Figure  4E). 

To  determine  if  stabilization  of  PI3K  contributes  to  the  oncogenic  potential  of  the 
mutants,  we  compared  overall  AKT  activity  to  total  pi  10a  levels  in  mutant  and  wildtype 
cells.  During  exponential  growth,  FI1047R  and  E545K-expressing  cells  had  3.9  and  2.4- 
fold,  respectively,  more  total  AKT  activity  than  those  expressing  wildtype  PIK3CA, 
though  the  corresponding  pi  1 0a  levels  were  only  2.0  and  1 .7-fold  higher  (Figure 
4A,F,G).  Given  that  these  mutant  forms  of  pi  10a  are  more  enzymatically  active  [15,  16], 
we  predict  that  the  combination  of  more  cells  in  the  pi  1 0ahl9h  state  and  higher  enzymatic 
activity  cooperate  to  induce  oncogenic  levels  of  AKT  activity. 

p110ahigh  cells  initiate  colony  formation 
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We  next  investigated  whether  pi  10ah'9h  cells  had  unique  functional  roles.  We  seeded 
single  cells  from  wildtype,  H1047R  and  E545K  parental  cell  lines  and  generated  clonal 
cell  lines.  Remarkably,  nearly  100%  of  the  clones  that  emerged  exhibited  a  bimodal 
profile  that  was  the  reverse  of  the  parental  population.  At  early  passage,  the  majority  of 
cells  were  in  the  pi  10ahigh  state  and  a  minority  in  the  pi  10alow  state,  with  corresponding 
changes  in  pAKT  (Figure  5A  and  Table  S2).  These  results  suggest  that  only  pi  1 0ah'9h 
cells  were  capable  of  forming  colonies,  at  least  partially  by  maintaining  high  AKT  activity. 
This  is  likely  important  for  overcoming  low  cell  density-related  and  replicative  stresses, 
as  we  observed  that  cell  lines  over-expressing  pi  10a,  either  wildtype  or  mutant,  survived 
many  more  serial  passages  than  clones  with  only  endogenous  levels  of  pi  10a  (Figure 
5B  and  Table  S2). 

If  enrichment  for  pi  1 0ah'9h  cells  is  important  in  the  early  stages  of  colony  formation,  we 
next  wanted  to  test  if  the  clones  reverted  back  to  the  parental  distribution  once  the 
population  established  steady  exponential  growth.  We  compared  pi  10a  levels  at  early 
and  late  passages  and  indeed  saw  reversion  to  the  parental  distribution  with  the  majority 
of  cells  in  the  pi  10alow  state  in  late  passages  (Figure  5C).  We  carefully  monitored  the 
rate  of  reversion  by  measuring  pi  10a  levels  in  sequential  passages  and  observed  a 
gradual  decrease  in  the  percent  of  cells  with  high  pi  10a  in  nearly  all  clones  by  passage 
8  (Figure  5D).  This  suggests  that  pi  1 0ahigh  cells  are  critical  for  colony  formation  and 
early  clonal  expansion,  however  heterogeneous  populations  of  predominantly  pi  10alow 
cells  facilitate  exponential  growth.  To  ensure  there  was  no  genetic  component  to  these 
effects,  we  established  a  second  generation  of  clones  from  primary  clonal  populations 
and  observed  the  same  trends  (Table  S2).  Additionally,  reverted  clonal  cell  lines 
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exhibited  similar  pi  10a  kinetics  after  EGF  stimulation  to  parental  populations  (Figure 
S4A-C). 

Interestingly,  at  passage  3  both  FI1047R  and  E545K-derived  clones  harbored  higher 
percentages  of  pi  10ahigh  cells  compared  to  wildtype  clones  and  reverted  to  the  parental 
distribution  at  a  slower  rate  than  wildtype  (Figure  5D).  This  is  consistent  with  our 
previous  observation  that  mutant  pi  1 0a  can  stabilize  the  pi  1 0ahigh  state  (Figure  4,  S4D). 
Importantly,  the  few  clones  that  failed  to  undergo  a  reversion  to  the  parental  distribution 
were  predominantly  mutant  clones  with  very  high  percentages  of  pi  1 0ahigh  cells  (Figure 
5D).  These  clones  ceased  to  proliferate  after  passage  3-4  and  exhibited  a  spindly 
morphology  and  senescence-associated  (3-galactosidase  activity  (Figure  5E).  We  find 
these  data  consistent  with  a  model  where  reversion  to  a  bimodal  distribution  of  primarily 
pi  10alow  cells  is  required  to  maintain  exponential  growth,  and  failure  to  revert  results  in 
oncogene-induced  senescence. 

Establishing  cell-cell  contacts  stabilizes  the  pi  10alow  state. 

If  low  cell  density  contributed  to  the  induction  of  high  pi  10a  during  colony  formation,  we 
wanted  to  explore  the  possibility  that  increasing  cell-cell  contacts  induces  the  pi  10alow 
state.  We  measured  pi  10a  levels  of  exponentially  growing  populations  every  24h  and 
found  that  pi  10a  and  pAKT  levels  steadily  decreased  as  cell  density  increased  (Figure 
6A,B).  The  decrease  in  total  pi  10a  can  be  attributed  specifically  to  a  drop  in  the 
proportion  pi  10ahigh  cells  and  an  accumulation  of  pi  10alow  cells  (Figure  6C).  This  is  also 
apparent  in  small  cell  clusters,  where  pi  10a  levels  are  heterogeneous  compared  to  lone 
cells  that  tend  to  be  pi  1 0ahigh  (Figure  6D).  Under  these  low-density  growth  conditions, 
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the  pi  10ahigh  population  is  enriched  for  cells  with  0-1  cell-cell  contacts,  while  the 
pi  1 0alow  population  is  enriched  for  cells  with  2-4  cell-cell  contacts  (Figure  6D).  Thus, 
cell-cell  contacts  may  modulate  pi  10a  levels  in  order  to  establish  contact 
inhibition/confluence  arrest. 

We  next  considered  the  possibility  that  cells  with  high  AKT  activity  are  enriched  at  the 
unencumbered  edges  of  tumors,  which  are  least  confluence-inhibited  and  highly 
proliferative.  These  “pushing  margins”  are  a  distinct,  though  molecularly 
uncharacterized,  feature  of  mutant  BRCAI-d riven  tumors  [25,  26].  We  thus  stained 
mammary  tumors  arising  from  MMTV-Cre;  BRCA1  flox/flox;  p53+/-  mice  for  pAKT  S473. 
We  detected  heterogeneous  pAKT-positivity  throughout  the  tumor  (Figure  S5A)  and  an 
enrichment  of  pAKT-positive  cells  along  pushing  margins  (Figure  S5B-F)  and  in  the 
lumens  of  ductal  hyperplasias  (Figure  S5G,FI).  Both  of  these  regions  represent  areas 
with  fewer  cell-cell  contacts,  supporting  our  hypothesis  that  low  cell  density  enhances 
the  pAKT-positive  state.  Flowever,  given  that  no  antibodies  against  pi  1 0a  or  p85a  are 
amenable  to  immunohistochemical  staining,  we  cannot  confirm  that  this  high  AKT 
activity  correlates  with  high  PI3K  protein. 
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Discussion 


Non-genetic  cell-to-cell  variability  has  been  observed  in  numerous  cellular  systems  and 
can  lead  to  distinct  cellular  fates.  The  most  classic  example  is  Ferrell  and  Machleder’s 
study  of  Xenopus  oocyte  maturation,  which  shows  that  variability  in  MAPK 
phosphorylation  leads  to  either  a  G2  or  metaphase  arrest  [27],  Similarly,  high  or  low 
levels  of  Nanog  determines  the  potential  of  embryonic  stem  cells  to  terminally 
differentiate  [10,  28],  and  high  or  low  levels  of  Sca-1  dictates  the  proclivity  for 
hemapotoietic  progenitor  cells  to  commit  to  the  erythroid  or  the  myeloid  lineage  [29], 
Another  cellular  fate  that  is  determined  by  variable  protein  expression  is  the  drug- 
tolerant/drug-sensitive  state.  In  certain  cell  lines,  gefitinib-resistance  can  be  conferred 
by  high  IGF-1R  signaling  and  high  KDM5A  expression  [7]  and  camptothecin-resistance 
by  high  DDX5  or  RFC1  expression  [4].  Our  study  now  shows  that  heterogeneity  in 
signal  transduction  is  another  meaningful  source  of  cell-to-cell  variability  that  results  not 
in  distinct  cell  fates,  but  rather  influences  the  transient  behavior  of  the  entire  population. 

We  show  that  a  bimodal  distribution  of  AKT  activation  is  an  invariable  characteristic  of 
exponentially  growing  MCF10A  cells.  We  propose  that  limiting  AKT  activity  to  only  20- 
30%  of  cells  in  a  population  serves  two  related  purposes:  1)  to  prevent  senescence;  and 
2)  to  maintain  sub-oncogenic  levels  of  PI3K  activity  in  large  populations. 

Our  data  show  that  clonal  populations  that  are  unable  to  revert  to  the  parental 
distribution  of  primarily  pi  10alow  cells  undergo  cellular  senescence.  In  these  clones, 
pi  10ah'9h;pAKTh'9h  cells  comprise  80-90%  of  the  population  immediately  prior  to 
senescence,  indicating  that  such  high  PI3K  activity  is  not  sustainable.  These  results  are 
consistent  with  models  of  p53-dependent  cellular  senescence  associated  with  the  loss  of 
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two  negative  regulators  of  PI3K,  PTEN  and  inositol  polyphosphate  4-phosphatase  type  II 
(INPP4B)  [30,  31].  Given  that  MCF10A  cells  are  p53-replete,  it  is  possible  that  over¬ 
accumulation  of  pi  1 0ahigh  cells  induces  senescence  by  the  same  mechanism. 

A  second  selective  advantage  of  maintaining  variability  in  PI3K  activity  may  be  to 
maintain  sub-oncogenic  levels  of  PI3K  activity  in  large  populations  such  as  tissues  and 
organs.  Given  that  over-activation  of  this  pathway  is  one  of  the  most  frequent  events  in 
cancer,  permitting  PI3K  activation  in  only  a  fraction  of  cells  within  a  population  may  be  a 
simple  mechanism  to  limit  overall  PI3K  activity.  However,  we  show  that  variability  is  not 
completely  overcome  in  populations  expressing  oncogenic  forms  of  pi  10a  or  in  epithelial 
tumor  sections,  which  may  have  important  implications  on  tumor  susceptibility  to  PI3K 
inhibitors. 

Assuming  that  only  tumor  cells  with  high  levels  of  pathway  activation  are  likely  to  die  in 
response  to  acute  inhibition  of  PI3K,  then  tumors  with  mosaic  patterns  of  pathway 
activation  (e.g.  as  judged  by  pAKT  staining)  will  exhibit  only  a  partial  response  to  a  PI3K 
(or  AKT)  inhibitor.  The  mosaic  pattern  of  pathway  activation  could  be  due  to  genetic 
variability  or  oscillation  between  high  and  low  levels  of  PI3K  expression  (as  seen  in 
MCF10A  cells).  In  the  latter  case,  it  may  be  possible  to  kill  a  larger  fraction  of  the  tumor 
cells  by  adjusting  the  therapy  to  ensure  that  the  subset  of  cells  that  are  protected  from 
PI3K  inhibitor  therapy  at  the  time  of  the  initial  dose,  become  exposed  to  the  inhibitor  a 
second  time  as  they  cycle  into  a  high  PI3K  pathway  state.  Thus,  a  detailed 
understanding  of  the  pharmacokinetics  of  the  drug  is  important  in  deciding  the  intensity 
and  frequency  of  dosing.  While  continuous  high  doses  of  PI3K  inhibitors  over  several 
days  or  weeks  might  effectively  kill  all  tumor  cells  as  they  cycle  into  the  high  activity 
state,  dose-limiting  toxicities  might  preclude  such  treatment.  An  alternative  approach  of 


26 


using  a  very  high  single  dose  once  a  week  could  avoid  toxicities  associated  with 
continuous  dosing  and  result  in  20-30%  tumor  cell  death  at  the  first  treatment,  followed 
by  a  period  in  which  another  20-30%  of  the  resistant  tumor  cells  cycle  into  the  high  PI3K 
pathway  state  and  die  at  the  second  treatment  and  etc. 

Various  pre-clinical  studies  using  PI3K  inhibitors  may  support  this  hypothesis.  Once 
daily  oral  administration  of  the  Novartis  pan-PI3K/mTOR  inhibitor,  NVP-BEZ235,  to 
tumor-bearing  mice  results  in  delayed  tumor  growth  but  rarely  tumor  regression.  pAKT 
is  ablated  in  these  tumors  1-2h  after  dosing,  but  reemerges  several  hours  later  [32-34], 
Similar  results  were  observed  using  the  Genentech  class  1 A  PI3K  inhibitor,  GDC-0941 
[35,  36].  Our  data  would  suggest  that  cells  that  were  PI3Klow;pAKTlow  during  the  initial 
dose  escaped  the  toxic  effects  of  the  drug,  and  cycled  to  the  PI3Khigh;pAKTh'9h  state 
several  hours  later  when  concentrations  of  the  drug  decreased.  The  possibility  that  PI3K 
inhibitors  are  selectively  targeting  the  advancing  edges  or  “pushing  margins”  of  tumors  is 
also  consistent  with  fact  that  these  tumors  are  not  advancing,  but  rather  stabilizing. 

As  we  emphasized  throughout  this  study,  variability  in  PI3K  activity  in  cell  populations  is 
a  regulated  process.  Though  the  mechanism  regulating  PI3K  degradation  and  re¬ 
synthesis  is  beyond  the  scope  of  this  study,  it  is  tempting  to  consider  inhibition  of  re¬ 
synthesis  to  dampen  PI3K  activity  in  all  cells.  Another  alternative  is  to  inhibit  the  E3 
ligase  that  targets  PI3K  for  degradation,  which  may  over-activate  the  pathway  in  all  cells 
and  induce  cellular  senescence.  This  “pro-senescence”  approach  has  recently  shown 
promise  in  a  prostate  cancer  xenograft  model  where  PTEN  is  pharmacologically 
inhibited  [37], 
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Materials  and  Methods 


Cell  lines  and  cell  culture 

MCF10A  cells  were  obtained  from  the  American  Type  Culture  Collection  and  cultured  as 
described  by  Debnath  et  a/ [38].  Stable  MCF10A  cell  lines  expressing  FIA-tagged 
wildtype  or  mutant  (FI1047R  or  E545K)  bovine  pi  10a  were  generated  by  retroviral 
infection  as  previously  described  [1 5].  Expression  of  wildtype  and  mutant  PIK3CA  cDNA 
was  driven  by  the  exogenous  CMV  promoter  of  the  JP1520  retroviral  vector  (J. 
Pearlberg).  Where  indicated,  cells  were  starved  in  growth  media  without  horse  serum, 
EGF  or  insulin  for  ~20h.  Acute  stimulation  with  growth  factors  was  done  by  adding 
growth  factor  directly  to  the  starvation  medium. 

Preparation  of  cells  for  flow  cytometry 

A  detailed  description  of  how  to  prepare  adherent  cells  for  analysis  by  flow  cytometry  is 
provided  in  the  Supplemental  Information.  Antibodies  used  for  flow  cytometry  were:  anti- 
HA.  1 1  -488  conjugate  (Covance);  anti-p-AKT  S473  (D9E)-488/647  conjugates,  anti-total 
AKT-647  conjugate,  anti-pEGFR  Y1068  (Cell  Signaling);  anti-p85a  (Upstate 
Biotechnology);  and  various  secondary  Alexa  Fluor  dyes  (Invitrogen).  Data  was 
analyzed  with  BD  FACSDiva  software  (BD  Biosciences),  FlowJo  software  (Tree  Star) 
and  Cytobank  [39]. 

Generation  of  clonal  populations 

Parental  pooled  populations  were  passaged  >9  times  in  selection  media  to  normalize  for 
infection  efficiency.  Parental  cells  were  trypsinized,  diluted  and  plated  at  a  density  of  1 
cell/well.  Colony  formation  was  monitored  daily  and  media  was  changed  every  three 
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days.  Once  colonies  grew  to  a  sufficient  size,  they  were  trypsinized  and  re-plated 
(passage  1)  and  serially  passaged  thereafter. 

Immunoassays 

Detailed  protocols  for  immunofluorescence,  immunoblotting,  and  immunohistochemistry 
are  provided  in  the  Supplemental  Information. 
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Figure  Legends 


Figure  1.  PI3K  pathway  activation  in  MCF10A  populations  is  heterogeneous. 

(A)  Cells  were  starved  and  acutely  stimulated  with  5nM  EGF  for  5’  with  or  without 
pretreatment  with  lOOnM  wortmannin.  Cells  were  then  harvested  and  stained 
with  anti-pAKT  S473  antibody  and  analyzed  by  flow  cytometry.  Unstimulated 
and  wortmannin  treated  cells  are  pAKT-negative  (P6).  EGF  stimulation  activates 
-30%  of  cells  and  shifts  them  into  the  pAKT-positive  gate  (P7).  Cells  stably 
expressing  the  oncogenic  FI1047R  pi  10a  mutant  exhibit  basal  AKT  activity  in  the 
unstimulated  state  (as  previously  reported),  and  EGF  stimulation  shifts  more 
cells  into  the  pAKT-positive  gate  (P7)  as  well  as  increases  the  amplitude  of  pAKT 
signal  (P8). 

(B)  Mammary  carcinomas  from  MMTV-Cre;  BFtCAIflox/flox;  p53+/-  mice  display 
heterogeneous  pAKT-positivity. 

(C)  MCF10A  cells  stably  expressing  FIA-tagged  pi  10a  were  treated  as  described  in 
(A)  and  co-stained  with  anti-pAKT  S473  and  anti-FIA  antibodies.  HA-pl  1 0a 
segregates  into  a  bimodal  distribution  and  correlates  with  pAKT  status  (Ft2=0.68). 

(D)  The  results  from  (C)  were  validated  by  immunofluorescence. 

Figure  2.  pi  10a  levels  are  dynamic  within  single  cells. 
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(A)  MCF10A  cells  stably  expressing  wildtype  HA-pl  10a  were  starved  and  acutely 
stimulated  with  5nM  EGF.  Cell  lysates  were  harvested  at  indicated  time  points 
and  analyzed  by  SDS-PAGE. 

(B)  MCF10A  cells  stably  expressing  wildtype  HA-pl  10a  were  treated  as  described  in 
(A).  Cells  were  harvested  at  indicated  time  points  and  analyzed  by  flow 
cytometry.  The  average  pi  10a  level  in  the  population  (solid  blue  line)  remains 
low  and  unchanging,  as  seen  by  western  blot  in  (A).  On  the  subpopulation  level, 
pAKT-positive  cells  (solid  red  line)  initially  have  high  levels  of  pi  10a,  which  drop 
precipitously  at  10’  and  recover  over  time.  The  pAKT-negative  subpopulation 
(dotted  red  line)  retains  constant  low  levels  of  pi  10a.  The  dynamic  changes  in 
pi  10a  levels  in  the  pAKT-positive  subpopulation  can  modestly  influence  pi  10a 
levels  of  the  entire  population,  as  seen  by  the  slight  increase  in  pi  1 0a  in  the 
presence  of  MG132  (dotted  blue  line). 

(C)  The  proportion  of  cells  in  the  pi  10ah'9h/med'um/low gates  was  measured  at  each  time 
point  in  a  representative  experiment  described  in  (B).  Fold  change  compared  to 
time  0  is  shown,  indicating  a  shift  of  cells  from  the  pi  10ahigh  to  pi  10amedlum  to 

pi  10alow  states  and  back.  The  short  4h  time  course  ensures  that  these  dynamic 
changes  are  not  due  to  cell  division. 

(D)  Average  AKT  activity  was  measured  on  the  population  and  subpopulation  levels 
in  the  cells  from  (B).  pi  10ahigh  cells  (solid  black  line)  maximally  activate  AKT 
compared  to  pi  10alow  cells  (dashed  black  line),  which  minimally  activate  AKT. 

pi  ioamedlum  cells  (dotted  black  line)  have  intermediate  levels  of  pAKT.  Given  that 
the  pi  10ahigh  cells  only  represent  -30%  of  the  population,  the  average  pAKT 
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level  on  the  population  level  (red  line)  is  only  slightly  higher  than  the  pi  10alow 
line. 


Figure  3.  Changes  in  pi  10a  levels  are  regulated  by  a  degradation/re-synthesis 
cycle. 

(A)  Asynchronously  growing  MCF10A  cells  expressing  wildtype  HA-pl  10a  were 
treated  for  4h  with  100|ig/mL  cyclohexamide  (CHX)  or  20|iM  MG  132  and 
analyzed  by  flow  cytometry.  (Left  panel)  On  the  population  level,  CHX  (blue  line) 
and  MG132  (orange  line)  treatment  shifted  the  pi  10a  distribution  to  the  left  and 
right,  respectively,  compared  to  untreated  cells  (green  line).  (Right  panel)  These 
shifts  in  pi  10a  levels  correlate  with  changes  in  overall  pAKT  level. 

(B)  On  the  subpopulation  level,  CHX  treatment  inhibits  the  accumulation  of  cells  in 
the  pi  1 0amedlum/hl9h  gates  by  retaining  cells  in  the  pi  1 0alow  state.  MG1 32 
treatment  results  in  the  accumulation  of  cells  in  the  pi  1 0amedlum/h'9h  gates. 
Changes  in  population  distribution  are  quantified  in  the  right  panel. 

(C)  Model  of  PI3K  degradation/re-synthesis  cycle  that  regulates  pathway  activation. 
(Left  panel)  Inactive  cells  are  pi  10alow;  pAKTlow  (Q4).  Extrinsic  signals  trigger  de 
novo  synthesis  of  PI3K  protein,  transitioning  cells  to  a  pi  10ah'9h;  pAKTlow  state 
(Q1),  where  they  are  competent  to  activate  the  pathway.  Once  a  growth  signal  is 
received,  cells  transition  to  the  pi  10ahigh;  pAKTh'9h  state  (Q2)  and  the  PI3K 
pathway  is  fully  activated.  Shortly  after  activation,  PI3K  is  rapidly  degraded  and 
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cells  enter  the  pi  10alow;  pAKThigh  state  (Q3).  Once  AKT  is  fully 
dephosphorylated,  cells  return  to  the  inactive  (Q4)  state.  (Right  panel) 
Immunofluorescence  of  EGF  stimulated  cells  shows  cells  in  each  of  these  four 
stages  of  the  cycle. 

Figure  4.  PIK3CA  mutations  stabilize  the  p110ah'9h  state. 

(A)  Asynchronously  growing  MCF10A  cells  stably  expressing  FIA-tagged  wildtype  or 
mutant  pi  10a  (FI1047R  or  E545K)  were  analyzed  by  flow  cytometry.  Cells 
expressing  the  FI1047R  or  E545K  mutations  had  2.0  and  1 .7-fold,  respectively, 
higher  levels  of  pi  10a  compared  to  wildtype. 

(B)  The  proportion  of  cells  within  the  pi  10ah'9h/medium/low  gates  was  measured  for  the 
three  cell  lines.  Both  mutant  cell  lines  had  more  cells  with  high  levels  of  pi  10a 
and  fewer  cells  with  low  levels  of  pi  10a. 

(C)  Asynchronously  growing  cells  expressing  wildtype  pi  10a  only  have  -25%  of 
cells  actively  signaling  or  primed  to  activate  signaling  compared  to  >40%  in 
either  mutant  cell  line.  Numbers  represent  an  average  of  three  experiments. 

(D)  Cells  expressing  wildtype  or  mutant  pi  10a  were  treated  as  described  in  Figure 
2B.  Total  levels  of  pi  1 0a  in  the  whole  population  did  not  change  over  time  (solid 
lines),  though  the  mutant  cells  retained  higher  average  pi  10a  levels  than  the 
wildtype  cells.  In  the  pAKTh'9h  subpopulations  (dashed  lines),  the  mutants 
followed  similar  dynamic  changes  in  pi  10a  levels  as  described  for  wildtype, 
though  with  much  greater  variability. 
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(E)  Asynchronously  growing  wildtype  and  mutant  lines  were  treated  with  CHX  or 
MG132  as  described  in  Figure  3A.  All  cell  lines  showed  a  decrease  in  pi  10a 
levels  following  treatment  with  CFIX.  Flowever  the  mutant  lines  were  less 
sensitive  to  the  MG132  treatment  than  the  wildtype  line. 

(F)  Average  AKT  activity  in  asynchronously  growing  wildtype  and  mutant  cell 
populations  was  measured  by  flow  cytometry.  The  FI1047R  and  E545K  mutant 
lines  exhibited  3.9  and  2.4-fold  increases,  respectively,  in  total  pAKT  level 
compared  to  wildtype. 

(G)  Wildtype  and  mutant  cells  were  analyzed  by  immunofluorescence  as  described 
in  Figure  1 . 


Figure  5.  Cells  with  high  levels  of  pi  10a  are  important  for  colony  formation. 

(A)  On  passage  3,  clonal  cell  lines  were  starved  and  acutely  stimulated  with  5nM 
EGF  and  analyzed  by  flow  cytometry.  In  each  set  of  histograms,  the  top  (blue) 
measures  pi  10a  and  the  bottom  measures  pAKT  S473  of  pi  10alow  (green)  and 
pi  10ah'9h  (red)  cells. 

(B)  Clones  over-expressing  wildtype  or  mutant  pi  10a  survived  more  passages  than 
cells  with  only  endogenous  levels  of  pi  10a.  Additionally,  expression  of  either 
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mutant  form  of  pi  10a  increased  clonal  viability  at  early  passages  compared  to 
expression  of  wildtype  pi  1 0a. 

(C)  Clones  at  late  passage  (>10)  were  analyzed  as  described  in  (A)  and  invariably 
reverted  to  the  parental  bimodal  distribution  characterized  by  a  large  pi  10alow 
subpopulation  and  a  small  pi  1 0ochigh  subpopulation. 

(D)  To  monitor  the  rate  of  reversion,  the  percent  of  cells  with  high  levels  of  pi  1 0a 
was  measured  by  flow  cytometry  for  six  serial  passages.  Mutant  clones  (middle 
and  bottom)  achieved  higher  pi  10ah'9h  percentages  and  reverted  slower  than 
wildtype  clones  (top).  Notably,  many  mutant  clones  (HR-5,6,7  and  EK- 
10,8,1 1 ,14)  failed  to  revert  and  instead  ceased  proliferating  at  early  passage. 

(E)  Clones  that  stopped  growing  stained  positively  for  senescence-associated  p- 
galactosidase  activity  compared  to  an  exponentially  growing  control  line. 

Figure  6.  Establishing  cell-cell  contacts  stabilizes  the  p110alow  state. 

(A)  Cells  achieve  ~1.5  doublings  every  24h,  thereby  increasing  the  number  of  cell¬ 
cell  contacts. 

(B)  Asynchronously  growing  cells  were  harvested  every  24h  and  analyzed  by  flow 
cytometry.  As  cell  density  increases,  average  pi  10a  (left)  and  pAKT  (right) 
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levels  in  the  population  decreases.  Data  shown  are  from  one  representative 
experiment. 

(C)  Analyzed  on  the  subpopulation  level,  the  decrease  in  pi  10a  can  be  attributed 
specifically  to  a  decrease  in  the  pi  10ahigh  subpopulation  and  an  accumulation  of 
the  pi  1 0alow  subpopulation. 

(D)  To  verify  these  results,  cells  were  grown  to  -50%  confluency  and  pi  10a  levels 
were  assessed  by  immunofluorescence  (Left  panel).  Cells  with  few  or  no  cell¬ 
cell  contacts  (white  circles)  stained  brightly  for  HA-pl  10a  whereas  cell  clusters 
with  many  cell-cell  contacts  (orange  boxes)  displayed  heterogeneity  in  HA-pl  10a 
staining  intensity.  The  number  of  cell-cell  contacts  was  counted  for  875 
individual  cells  and  compared  to  their  pi  1 0a  status  (Right  panel),  pi  1 0ah'9h  cells 
were  enriched  for  cells  with  0-1  cell-cell  contacts  and  pi  10alow  cells  were 
enriched  for  cells  with  2-4  cell-cell  contacts. 
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Figure  SI,  related  to  Figure  1 
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Figure  S2,  related  to  Figure  2 
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Mean  pi  10  level  (fold  change  over  Time  0) 


Figure  S3,  related  to  Figure  3 
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Figure  S4,  related  to  Figure  5 
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Figure  S5,  related  to  Figure  6. 
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Figure  SI,  related  to  Figure  1.  Correlation  among  other  PI3K  pathway  members. 


(A)  EGFR  activity  does  not  correlate  with  AKT  activity.  MCF10A  cells  were  serum  and 
growth  factor  starved  for  20h  and  acutely  stimulated  with  5nM  EGF.  (Left  panels)  Cells 
were  fixed  and  co-stained  with  anti-pAKT  S473  and  pEGFR  Y1068  antibodies  and 
analyzed  by  immunofluorescence.  Most  cells  display  some  degree  of  EGFR 
phosphorylation,  however  the  amplitude  of  pEGFR  staining  does  not  correlate  with 
pAKT  status.  (Right  panels)  Additionally,  cells  were  stimulated  with  EGF,  stained  with 
anti-pEGFR  Y1068  antibody,  and  analyzed  by  flow  cytometry.  EGF  stimulated  cells 
uniformly  shift  with  a  log-normal  distribution  to  the  positive  gate  (P9),  in  contrast  to  the 
bimodal  distribution  observed  for  pAKT  S473. 

(B)  Total  AKT  level  does  not  determine  AKT  activity.  MCF10A  cells  were  treated  as 
described  in  (A).  Cells  were  harvested  and  co-stained  with  anti-pAKT  S473  and  anti¬ 
total  AKT  antibodies  and  analyzed  by  flow  cytometry.  (Top  right  panel)  Total  AKT  levels 
are  log-normal  in  MCF10A  populations,  compared  to  pAKT  S473  levels  (left  panel), 
which  are  bimodal.  pAKT-positive  cells  (P8)  have  slightly  higher  levels  of  total  AKT  than 
pAKT-negative  cells  (P7),  however  there  is  significant  overlap  between  the  two 
populations,  weakening  the  correlation  (R2=0.59). 

(C)  Bimodal  AKT  activation  is  not  ligand-specific.  MCF10A  cells  were  serum  and  growth 
factor  starved  and  acutely  stimulated  with  2nM  EGF  or  1  |iM  insulin  for  10  minutes.  Cells 
were  harvested  and  stained  for  pAKT  S473  and  analyzed  by  flow  cytometry.  Similar 
bimodal  distributions  of  AKT  activity  were  obtained  with  EGF  and  insulin  stimulation, 
suggesting  that  heterogeneous  AKT  activation  is  not  a  ligand-specific  effect. 
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(D)  p85a  level  correlates  with  pi  1 0a  level.  MCF1 OA  cells  expressing  HA-pl  1 0a  were 
harvested  and  co-stained  with  anti-HA  and  anti-p85a  antibodies.  (Top  right  panel) 
Endogenous  p85a  levels  segregate  into  a  bimodal  distribution,  similar  to  HA-pl  10a 
levels  (left  panel),  pi  10ahigh  cells  (P6)  are  also  p85ah'9h  and  pi  10alow  cells  (P5)  are 
p85alow,  conferring  a  strong  correlation  (R2=0.71). 

(E)  Two  non-transformed  mammary  epithelial  cell  lines,  MCF10A  and  HMEC,  and  three 
transformed  breast  cancer  cell  lines,  SUM159PT,  HS578T  and  MDA-MB-231 ,  were 
serum  starved  for  24  hours  and  acutely  stimulated  with  saturating  doses  of  growth 
factors  (5nM  EGF  +  luM  insulin)  for  10  minutes  to  maximally  activate  endogenous  PI3K 
within  the  cell  population.  Cells  were  harvested  as  described  in  the  Methods  and 
analyzed  by  flow  cytometry  for  phosphorylation  of  AKT  at  S473.  Dotted  lines  represent 
AKT  activity  in  starved  cells  and  solid  lines  represent  AKT  activity  in  stimulated  cells. 
Only  non-transformed  cells  exhibited  distinct  bimodality,  while  the  transformed  cells 
exhibited  a  more  complete  shift  to  AKT  positivity.  This  suggests  that  in  non-transformed 
cells  AKT  activity  is  negatively  regulated,  perhaps  via  the  modulation  of  PI3K  protein 
levels  as  described  in  this  manuscript.  In  order  to  achieve  oncogenic  levels  of  AKT 
activity,  cancer  cells  may  override  this  negative  regulation  by  maintaining  constitutively 
high  levels  of  PI3K  or  perhaps  by  some  means  of  lowering  the  threshold  for  PI3K-AKT 
activation. 

Figure  S2,  related  to  Figure  2.  p85a  kinetics  closely  mirror  pi  10a  kinetics  during  a 
timecourse  of  EGF  stimulation. 

(A)  MCF10A  cells  were  serum  starved  and  stimulated  with  5nM  EGF.  p85a  and  HA-pl  10a 
protein  levels  were  measured  at  various  time  points  by  flow  cytometry  using  antibodies 
against  endogenous  p85a  and  HA.  As  shown  in  Figure  2A,  HA-pl  10a  levels  in  the 
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pAKT-high  population  drop  precipitously  at  10  minutes.  This  drop  is  mirrored  by 
endogenous  p85a  levels,  which  drop  to  a  similar  but  slightly  lesser  degree.  This  may  be 
explained  by  the  fact  that  there  is  more  p85a  than  pi  10a  in  the  cell,  and  monomeric 
p85a  may  not  be  subject  to  the  same  degradation  mechanism  as  heterodimerized  p85a. 

(B)  The  distribution  of  cells  within  the  p85alow,  p85amedlum  and  p85ah'9h  populations  closely 
resembles  the  distribution  within  pi  1 0a  populations,  including  the  drop  in  the  proportion 
of  p85ahigh  cells  at  10  minutes,  which  correlates  with  the  drop  in  p85a  protein  levels  in 
(A). 

(C)  The  AKT  activity  of  cells  within  the  p85alow,  p85amedlum  and  p85ah'9h  populations  closely 
resembles  those  of  pi  1 0a  populations,  with  the  p85ahigh  population  driving  maximal  AKT 
activity  within  the  total  population. 

Figure  S3,  related  to  Figure  3.  MG132  treatment  can  stabilize  pi  10a  levels  during  EGF 
stimulation. 

(A)  Cells  were  serum  starved  for  24  hours  and  pre-treated  with  50uM  MG132  for  1h  prior  to 
5nM  EGF  stimulation,  pi  1 0a  levels  were  measured  at  various  time  points  by  flow 
cytometry.  MG  132  treatment  resulted  in  a  partial  rescue  of  pi  10a  degradation  at  10 
minutes  and  an  overall  increase  in  pi  10a  protein  in  the  total  cell  population.  We  were 
unable  to  achieve  a  complete  rescue  of  pi  10a  degradation,  which  may  be  due  to  the  fact 
that  we  were  unable  to  treat  cells  with  higher  concentrations  of  MG132  or  for  longer 
durations  without  greatly  altering  AKT  activity.  It  is  also  possible  that  pi  10  degradation 
is  carried  out  by  cellular  proteases  in  addition  to  the  proteosome. 
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(B)  50uM  MG132  treatment  was  also  sufficient  to  block  the  drop  in  the  proportion  of  cells  in 
the  pi  10ahigh  state  at  10  minutes  that  is  observed  in  untreated  cells. 

(C)  The  raw  data  from  the  24h  timepoint  in  the  experiment  described  in  Table  SI  are  plotted 
as  dot  plots  to  demonstrate  little  to  no  correlation  between  any  parameter  parings 
besides  pAKT  and  pi  10a  protein. 

Figure  S4,  related  to  Figure  5.  Late  passage  clonal  cell  lines  behave  similarly  to  parental 
pooled  cell  lines  after  reversion  to  the  parental  bimodal  distribution. 

(A)  Two  late  passage  clonal  cell  lines  and  one  parental  cell  line  were  serum  starved  for  24 
hours  and  acutely  stimulated  with  5nM  EGF.  pi  1 0a  levels  were  measured  at  various 
time  points  by  flow  cytometry.  Despite  having  the  reversed  bimodal  profile  at  early 
passage,  late  passage  clonal  cell  lines  that  have  reverted  to  the  parental  distribution 
exhibit  the  same  pi  1 0a  kinetics  as  parental  cell  lines. 

(B)  Population  distributions  in  late  stage  clonal  cell  lines  resemble  the  population  distribution 
of  a  parental  cell  line,  including  the  drop  in  the  proportion  of  cells  in  the  pi  1 0ah'9h  state  at 
10  minutes. 

(C)  In  clonal  cell  lines,  the  AKT  activity  of  cells  within  the  pi  10alow,  pi  10amedlum  and  pi  10ah'9h 
populations  closely  resembles  those  of  a  parental  cell  line,  with  the  pi  1 0ahigh  population 
driving  maximal  AKT  activity  within  the  total  population. 

(D)  Clonal  cell  lines  derived  from  WT,  H1047R  or  E545K  parental  pooled  cell  lines  were 
harvested  at  intermediate  passages  and  analyzed  for  pi  10a  levels  and  AKT  activity 
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following  serum  starvation  and  acute  5nM  EGF  stimulation.  As  we  showed  for  parental 
cell  lines  in  Figure  4A,  clonal  cells  lines  expressing  mutant  PIK3CA  have  higher  levels  of 
pi  10a  and  pAKT  than  clonal  cell  lines  expressing  wildtype  PIK3CA. 

Figure  S5,  related  to  Figure  6.  pAKT  is  enriched  at  pushing  margins  and  ductal 
hyperplasias. 

(A)  Mammary  tumors  from  MMTV-Cre;  BRCA1  flox;flox;  p53+/-  mice  were  stained  for  pAKT 
S473.  Near  the  center  of  the  tumor,  pAKT-positive  cells  are  heterogeneously 
distributed,  as  indicated  by  arrowheads. 

(B-F)  pAKT  is  enriched  at  the  non-invasive,  smooth  edges  of  these  tumors,  referred  to  as 

“pushing  margins”.  Arrowheads  roughly  enclose  these  margins  in  five  separate  MMTV-Cre; 

BRCA1  flox;flox;  p53+/-  tumors. 

(G-H)  pAKT  is  also  enriched  in  the  lumens  of  ductal  hyperplasias  (arrowheads). 
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Table  SI,  related  to  Figure  3. 


R-squared  value 

Parameter  pairing 

24h 

48h 

72h 

96h 

pi  1 0,  pAKT 

0.523 

0.536 

0.596 

0.447 

pi  1 0,  DNA  content 

0.017 

0.004 

0.003 

0.006 

pi  10,  Cell  size 

0.030 

0.156 

0.155 

0.171 

pi  10,  Cell  granularity 

0.332 

0.277 

0.280 

0.330 

pAKT,  DNA  content 

0.033 

0.091 

0.125 

0.260 

pAKT,  Cell  size 

0.144 

0.046 

0.038 

0.005 

pAKT,  Cell  granularity 

0.262 

0.193 

0.220 

0.209 

DNA  content,  Cell  size 

0.127 

0.172 

0.194 

0.206 

DNA  content,  Cell  granularity 

0.044 

0.046 

0.032 

0.026 

Cell  size,  Cell  granularity 

0.357 

0.384 

0.369 

0.350 
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Table  S2,  related  to  Figure  5. 


Parental  Population 

Colonies  formed 

Survival  to  P2 

P4  Clones  with  high  pi  10a 

MCF10A 

20/48  (0.42) 

4/20  (.20) 

ND 

PIK3CA  WT 

24/48  (0.50) 

16/24  (0.67) 

10/11  (.91) 

PIK3CA  H1047R 

22/64  (0.34) 

18/22  (0.82) 

16/16  (1.0) 

PIK3CA  E545K 

27/72  (0.38) 

24/27  (0.89) 

22/22  (1 .0) 

Clone  HR-2 

9/24  (.38) 

6/9  (.67) 

5/6  (.83) 

Clone  EK-3 

9/24  (.38) 

7/9  (.78) 

1/6  (.17)* 

ND  -  Not  determined;  *  Cells  senesced  before  enough  could  be  harvested  for  analysis 
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Table  SI,  related  to  Figure  3.  Correlation  coefficients  for  various  parameter  pairings. 


Asynchronously  growing  E545K  cells  were  harvested  24h,  48h,  72h,  and  96h  after  initial  plating. 
Cells  were  stained  with  anti-pAKT  S473  and  anti-HA  antibodies  and  propidium  iodide  (PI)  to 
monitor  DNA  content.  Cells  were  analyzed  by  flow  cytometry  and  the  raw  data  from  >10,000 
cells/sample  were  collected  for  the  following  parameters:  pi  10a  protein,  pAKT  status,  DNA 
content  (PI),  cell  size  (FSC),  cell  granularity  (SSC).  The  correlation  coefficients  (R2  value)  for 
the  indicated  parameter  pairings  were  calculated,  demonstrating  that  pi  10a  protein  level-pAKT 
status  is  the  only  parameter  pair  with  a  strong  correlation. 

Table  S2,  related  to  Figure  5.  Statistics  from  the  generation  of  clonal  populations.  Single 
cell  clones  were  derived  from  uninfected  MCF1 0A,  pooled  wildtype  pi  1 Oa-expressing,  pooled 
mutant  pi  10a  -expressing,  and  mutant  clonal  parental  populations. 
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Supplemental  Methods 


Preparation  of  cells  for  flow  cytometry 

Cells  were  treated  as  described  in  Figure  Legends.  Cells  were  then  washed  once  with  PBS  and 
gently  lifted  from  the  plate  using  a  cell  scraper  with  a  flexible  rubber-like  tip  (Sarstedt)  in  1  mL 
PBS.  Cells  were  briefly  triturated  with  a  pipette  to  obtain  a  single-cell  suspension  and 
immediately  fixed  by  adding  PFA  to  a  final  concentration  of  3.7%.  Cells  were  rocked  at  room 
temperature  for  10  minutes,  spun  down  and  transferred  to  5mL  polypropylene  tubes.  Cells  were 
then  permeabilized  by  adding  3mL  100%  cold  methanol  drop-wise  while  vortexing  at  medium 
speed.  Tubes  were  then  transferred  to  ice  for  >1  h  or  stored  at  -20 °C  for  up  to  two  weeks  before 
antibody  staining. 

For  staining,  cells  were  washed  twice  with  3mL  staining  buffer  (0.5%  BSA  +  0.02%  sodium 
azide  in  PBS)  and  rehydrated  by  10  minutes  incubation  at  room  temperature  in  staining  buffer 
after  each  wash.  Approximately  1x10s  cells  were  stained  per  sample  in  a  total  volume  of  IOOuL 
staining  buffer  plus  antibody  cocktail.  Cells  were  labeled  for  30-60  minutes  at  room  temperature 
in  the  dark.  Cells  were  washed  twice  in  staining  buffer  and  stained  with  secondary  antibodies  if 
necessary.  Cells  were  then  resuspended  in  PBS  and  analyzed  on  a  BD  LSR  II  instrument  (BD 
Biosciences).  Data  was  analyzed  with  BD  FACSDiva  software  (BD  Biosciences),  FlowJo 
software  (Tree  Star)  and  Cytobank.  Gates  are  delineated  by  negative  controls  (no  EGF  or 
wortmannin  treatment)  in  the  case  of  pAKT  plots.  In  H A-pl  10  plots,  gates  are  set  in  regions 
where  it  is  apparent  that  there  are  distinct  modes. 

Immunofluorescence 

Cells  were  grown  on  glass  coverslips  and  treated  as  described  in  figure  legends.  Cells  were 
fixed  in  3.7%  PFA,  washed  with  PBS  and  permeabilized  in  100%  cold  methanol  at  -20 °C.  Cells 
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were  washed  and  blocked  with  1%  horse  serum  in  PBS  for  30  minutes  at  room  temperature. 
Primary  antibody  staining  was  done  for  1  h  at  room  temperature  followed  by  two  PBS  washes 
and  secondary  antibody  staining  for  1h  at  room  temperature  in  the  dark.  Antibodies  used  for 
immunofluorescence  were:  anti-pAKT  S473  (193H12),  anti-pEGFR  Y1068  (Cell  Signaling);  anti- 
HA.1 1  (Covance);  and  various  secondary  Alexa  Fluor  dyes  (Invitrogen). 

Immunoblotting 

Cells  were  lysed  in  RIPA  buffer  (50mM  Tris  pH  7.4,  150mM  NaCI,  2mM  EDTA,  1%  NP-40,  0.1% 
SDS,  0.4%  sodium  deoxycholate)  supplemented  with  leupeptin  (4pg/mL),  aprotinin  (4pg/mL), 
pepstatin  (4pg/mL),  PMSF  (ImM),  NaF  (50m M),  NaPPi  (5mM),  p-glycerophosphate  (10mM), 
Na3V04  (5mM)  and  DTT  (0.5mM).  Lysates  were  cleared  by  centrifugation  at  1 4,000  x  g  for  1 0 
minutes  at  4°C,  and  the  supernatants  were  subjected  to  SDS-PAGE  and  transferred  to 
polyvinylidene  difluoride  membranes.  Immunoblotting  was  performed  with  the  following 
antibodies:  anti-FIA.11  (Covance),  anti-pi  10a  (BD  Transduction  Laboratories),  anti-p85a 
(Upstate  Biotechnology),  anti-pAKT  S473  (Cell  Signaling)  and  anti-p-actin  (Abeam). 

Immunohistochemistry 

MMTV-Cre;  BRCA1  flox/flox;  p53+/-  mice  were  bred  and  maintained  in  accordance  with  Beth 
Israel  Deaconess  Medical  Center  IACUC  guidelines.  Tumors  were  fixed  in  SafeFix  II  solution 
(Fisher  Scientific)  and  paraffin  embedded.  Immunohistochemical  staining  with  anti-pAKT  S473 
(IHC  specific)  (Cell  Signaling)  was  performed  according  to  the  antibody  manufacturer’s 
recommendations  with  the  following  exceptions:  antigen  unmasking  was  done  in  citrate  buffer  at 
a  sub-boiling  temperature  for  20  minutes;  and  PBST  was  used  as  the  wash  buffer. 
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Appendix  B 


Bibliography  and  Personnel  List 
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